Shortcut to a Fermi-Degenerate Gas of Molecules via Cooperative Association 
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We theoretically examine the creation of a Fermi-degenerate gas of molecules by considering a 
photoassociation or Feshbach resonance applied to a degenerate Bose-Fermi mixture of atoms. This 
problem raises an interest because, unlike bosons, fermions in general do not behave cooperatively, 
so that the collective conversion of, say, two million atoms into one million molecules is not to 
be expected. Nevertheless, we find that the coupled Fermi system displays collective Rabi-like 
oscillations and adiabatic passage between atoms and molecules, thereby mimicking Bose-Einstein 
statistics. Cooperative association of a degenerate mixture of Bose and Fermi gases could therefore 
serve as a shortcut to a degenerate gas of Fermi molecules. 
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Laser cooling of molecules is difficult, if not impos- 
sible, to apply in practice. This bottleneck occurs be- 
cause the removal of energy from the translational de- 
grees of freedom is most often accompanied by excita- 
tion of ro-vibrational modes which, upon de-excitation, 
couple back into the translational motion in a manner 
that heats the system. A laser-cooling-based approach 
[1] to a Fermi-degenerate gas of molecules is therefore 
unlikely to succeed. Moreover, buffer-gas cooling [2] and 
Stark-deceleration [3] techniques, while no doubt of prac- 
tical interest in their own right, currently view quantum 
degeneracy from the horizon. 

One possible route to ultracold molecules is photoas- 
sociation [4], and the formation of nondegenerate Fermi 
molecules was indeed recently observed [5]. A shortcut 
to quantum degenerate molecules is in the early stages 
of development for Bose-Einstein condensates (BECs) , so 
far delivering strongly enhanced molecule formation [6], 
precise measurements of the light-induced shift of the 
binding energy [7], and searches for a fundamental limit 
to the atom-molecule conversion rate [8,9]. The mag- 
netic field counterpart to photoassociation is the Fesh- 
bach resonance [10]. Often seen as a tool for adjusting 
the scattering lengths, thereby enabling condensation in 
otherwise-incondensible systems [11,12], the latest stud- 
ies of a BEC tuned near a Feshbach resonance have cul- 
minated in evidence for atom-molecule coherence [13]. 

Motivated by these exciting experiments, we report 
that a degenerate gas of Fermi molecules could be formed 
by applying either a photoassociation [5-9] or Feshbach 
resonance [11-13] to a degenerate mixture of Bose and 
Fermi atoms [14-16]. The collective Fermi atom-molecule 
system is predicted to undergo Rabi-like oscillations, 
as well as adiabatic passage, from atoms to molecules, 
thereby mimicking cooperative behavior that was previ- 
ously attributed to Bose statistics [17,18]. 

Photoassociation occurs when two atoms absorb a laser 
photon, thereby jumping from the free two-atom contin- 
uum to a bound molecular state. If the initial atoms 
form a Bose-Einstein condensate, then a photoassocia- 



tion laser could conceivably be used to convert a BEC of 
atoms into a BEC of molecules [19-21]. Collective free- 
bound photoassociation is theoretically identical to mag- 
netoassociation, whereby a molecular condensate can be 
created when one atom from a condensed pair spin flips 
in the presence of a magnetic field tuned near a Feshbach 
resonance [22-24]. Intuition developed in one instance is 
therefore generally applicable to the other, and we will 
often refer simply to collective association. 

Analogous with coherent optical transients in few level 
atomic systems [25] , photoassociation of a BEC has been 
predicted to induce Rabi-like oscillations between atomic 
and molecular condensates [19-21], whereby an entire 
gas of, say, two million Bose-condensed atoms are col- 
lectively converted into a million molecules that are, in 
turn, collectively converted back into (roughly) two mil- 
lion atoms, ad infinitum. Another interesting possibility 
arises because the ground state of the system is all atoms 
for large positive detunings (far below threshold) and all 
molecules for large negative detunings (far above thresh- 
old), so that a slow sweep of the laser detuning from one 
extreme to the other will collectively convert a BEC of 
atoms into a BEC of molecules [19]. Incidentally, it was 
a particular combination of these two concepts, applied 
instead to magnetoassociation, that led to the observa- 
tion [13] of collective Ramsey fringes between an atomic 
condensate and a small fraction of molecular condensate 
dressed by dissociated atom pairs [26-28] . 

The statistics of neutral atoms is determined by the 
number of neutrons in the nucleus, which must be odd 
for fermionic atoms and even for bosonic atoms. The 
sum of the total number of neutrons in the nucleus of 
the constituent atoms similarly determines the statistics 
of neutral molecules. Accordingly, molecules formed by 
free-bound association of two fermions will necessarily 
result in a boson, whereas fermionic molecules are born 
from the union of a boson and a fermion. Given a de- 
generate mixture of Bose and Fermi gases [14-16], is it 
possible that collective free-bound association could serve 
as a source of degenerate Fermi molecules? 
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To address this question, we model a degenerate 
Bose-Fermi mixture of atoms [14-16] coupled by ei- 
ther a Feshbach or photoassociation resonance to a 
Fermi-degenerate gas of molecules. The initial bosonic 
[fermionic] atoms are denoted by the field <f>(r,t) 
[ip-(r, t)], and the fermionic molecules by the field 
il>+(r,t). Neglecting particle-particle interactions, the 
Hamiltonian density for such an untrapped system can 
be written H = Ho + Hi, where 
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Here ms (to_) is the mass of a bosonic (fermionic) atom, 
m + = ms + m— is the mass of a molecule, Hb{±) is the 
so-called chemical potential that implicitly accounts for 
particle trapping, So is the binding energy (detuning) of 
the molecular state relative to the dissociation threshold, 
and JC is the (real) atom-molecule coupling. 

We now make some simplifications to allow for ease of 
modeling. First, we scale the Fermi (Bose) atom fields 



e lfJ - Bt 6), and the molecule 



as tjj- — > e~ 

field as ip+ — > e^ l ^- +IJ - B ^ t ip + , which yields /i+ — > fi' + = 
/i+ — (/i_ + /Lis); in turn, is absorbed into the detun- 
ing. Second is that atom-molecule conversion occurs on a 
timescale much faster than the motion of the atoms in the 
trap, allowing us to neglect the kinetic energies, and jus- 
tifying our omission of an explicit trapping potential for 
the particles. Third is that, whatever the Fermi energy 
associated with the fermionic components, it lies within 
the Wigner threshold regime, so that the coupling K. can 
be taken as the same for all modes. Finally, we switch 
to momentum space, but retain only the k = atomic 
condensate modes since, due to Bose stimulation, these 
transitions are favored over k ^ modes. The simplified 
Hamiltonian reads 



H = \$ c L c k _ ^K(clb k a + a^Ck) 
k 



(2) 



where a annihilates a bosonic atom with wavevector 
k = 0, 6k (ck) annihilates an atom (molecule) with 
wavevector k, and k = JC/VV. 

The system described by the Hamiltonian (2) is a par- 
ticular nonlinear version of a two-level atom driven by an 
external field. Of course, the two- level atom is well known 
to undergo Rabi flopping [25], and similar behavior ob- 
served for the present system would in principle answer 
question of whether collective association will provide a 
shortcut to the formation of a Fermi-degenerate gas of 
molecules. We therefore consider an on-resonance sys- 
tem, for which S = and the Heisenberg equations of 
motion are 
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Z&k = KCl^Ck, 

icu = nb^a. 
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If the system starts out with Nb bosonic atoms, the 
operator a will have a characteristic size \/Nb. On the 
other hand, while the fermionic operators 6k and Ck are 
characteristically of order unity, the sum over momen- 
tum states will ultimately introduce the initial number of 
fermionic atoms. Hence, for equal numbers of bosons and 
fermions, Nb = Np = N, it is intuitively obvious that 
the equations of motion (3) yield a system that evolves 
with the characteristic frequency Q — VN k. 

Using Fock states, said intuition is confirmed by solv- 
ing the time-dependent Schrodinger equation numeri- 
cally, the results of which are shown in Fig. 1. For the 
simplest case of N = 1, the system undergoes complete 
conversion from a doubly-degenerate Fermi-Bose system 
of atoms to a Fermi-degenerate gas of molecules in a time 
Q^ 1 . However, considering a larger number of initial 
atoms, N — 5, we see that quantum many-body fluctu- 
ations not only frustrate complete conversion, but also 
adjust the oscillation frequency and lead to collapse and 
revival. Increasing the initial particle number to TV = 10 
brings the amplitude of the first half-oscillation closer 
to the ideal case. This behavior is exactly analogous to 
the single-component bosonic case [19,21]. Although lim- 
ited computational resources preclude the explicit inves- 
tigation of large particle number, based on the bosonic 
analogy [19,21] we fully expect the first half oscillation 
to be complete for large N, i.e., N ^> 1 initial Fermi 
atoms should be converted to N Fermi molecules over a 
timescale r <~ il^ 1 . 

Nevertheless, collisions between particles will shift the 
energy of a state with a given number of atoms and 
molecules in a manner that depends nonlinearly on the 
numbers, making it difficult to maintain exact resonance. 
Furthermore, once a the system has Rabi flopped to 
molecules, the depopulated atomic condensate makes 
Bose stimulation, and thus neglect of the k / BEC 
modes, questionable. And so we investigate the more ro- 
bust possibility of rapid adiabatic passage [25] . From the 
Hamiltonian (2), it should be clear that the system will 
favor all atoms for large positive detunings, while favor- 
ing all molecules for large negative detunings. With fi 
established as the characteristic frequency for collective 
atom-molecule conversion, changes in the detuning that 
are of the order of £1, and occur over a time of order 
should qualify as adiabatic. Hence, if the detun- 
ing is swept in a linear fashion according to S = —££l 2 t, 
then dimensionless sweep rates £ ~ 1 should enable rapid 
adiabatic passage to a degenerate gas of Fermi molecules. 

Our suspicions are again corroborated by a Fock-state- 
basecl numerical solution to the Schrodinger equation, 
shown in Fig. 2. While many-body effects appear to 
rather seriously affect the efficiency of a marginally adi- 
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abatic sweep (£ = 1) compared to the N = 1 case, the 
difference between TV = 5 and TV = 10 (not shown) is in 
fact small. Overall, many-body effects are expected to 
be weak for near-adiabatic sweeps (£ ~ 1), and vanish- 
ing for sweeps that are deep-adiabatic (£ ~ 0.1), again in 
agreement with our BEC experience [19]. On the mat- 
ter of photodissociation to noncondensate modes-and the 
related pair correlations [26-28]-we note that, in the all- 
boson case [29], such transitions can be neglected for a 
sweep directed as in Fig. 2, i.e., for 5 < 0. Given the 
success of the analogy so far, similar results are expected 
for a Fermi system. 

Before closing, we estimate some explicit numbers. 
We eschew photoassociation because of the losses associ- 
ated with the electronically-excited state [30], and focus 
on the atom-molecule coupling provided by the 87 Rb- 
40 K Feshbach resonance located at Bq = 534 G [31], 
which has a width Ar = 4 G and an associated 
zero-field Fermi-Bose-atom scattering length a^a = 
— 17.8nm. Accordingly, the atom-molecule coupling is 
K « (47r|af f f|^BohrAi?/m r ) 1 / 2 = 0.14cm 3 / 2 s" 1 , where 
we have estimated the difference between the Fermi-Bosc 
atom pair and molecular magnetic moments to be equal 
to the Bohr magneton /XBohr- Assuming Nb — 10 5 con- 
densate atoms in a trap with respective radial and axial 
frequencies u r /2ir = 215 Hz and oj a /2n — 16.3 Hz [16], 
the density of bosons is ps = 8.1 x 10 13 cm~ 3 . As for the 
fcrmions, we assume a modest number, say, Np = 10 , 
which has three consequences: (i) the atomic BEC will 
act as a reservoir, thus absorbing any heat created by 
holes in the Fermi sea [32]; (ii) barring an unfortu- 
nately large scattering length for Bose-atom and Fermi- 
molecule collisions, we can neglect the possibility of any 
Fermi-Bose collapse instabilities [16]; (hi) the size of the 
Fermi cloud (Rp — 8.3 /im) is smaller than the BEC 
(Rb = 10 //m), so that overlap is not an issue. Moreover, 
for Nb 3> Np, the timescale for atom- molecule conver- 
sion is T a 2 m ~ (y/pij/C) -1 = 8.2 x 10~ 7 s. This timescale 
is safely below the fastest timescale for trapped-atom mo- 
tion r t = (u) r /2ir)~ 1 = 4.7 x 10 -3 s, justifying our neglect 
of trap dynamics and the kinetic energy; physically put, 
this means that the Fermi energy is negligible compared 
to the atom-molecule coupling strength. 

Lastly we discuss the neglected role of particle-particle 
interactions, i.e., collisions. These are described in terms 
of a single parameter, the s-wave scattering length a, 
where the value of a differs depending on whether one has 
in mind Bose-Bose atom, Fermi-Bose atom, Fermi-Bosc 
atom-molecule, or Fermi-Fermi atom-molecule collisions 
(s-wave Fermi-Fermi atom and molecule collisions are of 
course Pauli blocked). Whatever the species involved, 
the real part of a describes elastic collisions between par- 
ticles, and leads to so-called mean-field shifts. The imagi- 
nary part need only be considered for atom-molecule col- 
lisions terms, and arises because the molecules are vi- 
brationally very hot, and therefore undergo vibrational 
relaxation induced by collisions with unconverted atoms. 



The collisional coupling strength is then A = 2tt ha/m*, 
where m* is the mass of the atom or the reduced mass of 
the atom-atom (atom-molecule) pair, and it is included in 
the theory as a density dependent detuning. The units 
of the collisional coupling are [A] = cm 3 /s, so that an 
appropriate density factor is needed to arrive at a col- 
lisional rate. Estimating the real part of a from the 
87 Rb- 40 K collisions [31], and the imaginary part from 
collisions [22,24], a typical density (p ~ 10 14 ) leads to 
elastic and inelastic collisions that occur on a timescale 
t c <~ (/oA) _1 ~ 1ms. For rapid adiabatic passage, the 
contribution from the collisional interaction is negligible 
compared to the detuning, except near-resonance; how- 
ever, the system only spends about T a i m ~ in this 
region, which is short enough to expect collisions to be 
negligible. If, by chance, collisions do turn out to be a 
problem, the trap could be expanded to further reduce 
the rate of collisions with unconverted atoms. 

In conclusion, we highlight that the term Q = >/N k 
was previously referred to as the Bose-enhanced free- 
bound coupling, and the detuning sweep, for example, 
referred to as Bose-stimulated rapid adiabatic passage 
from atoms to molecules. However, this behavior is now 
clearly independent of statistics, so that Bose stimulation 
of free-bound association has nothing whatsoever to do 
with Bose statistics, but is instead a many-body cooper- 
ative effect that applies equally well to Fermi-degenerate 
systems. We expect this situation to arise whenever 
the system is addressed as a unit, i.e., when the atom- 
molecule coupling strength is larger than the Fermi en- 
ergy. To four-wave mixing in a Fermi gas of atoms [33,34], 
we therefore add a further example of fcrmions mimicking 
an effect that was previously attributed to Bosc-Einstein 
statistics. Such behavior should allow for a shortcut to 
a degenerate gas of Fermi molecules via a collective pho- 
toassociation or Feshbach resonance. 

The authors acknowledge Juha Javanaincn for helpful 
discussions. This work supported by the Academy of 
Finland (MM and KAS, project 50314), and the Magnus 
Ehrnrooth foundation (OD). 



[1] B. DeMarco and D. S. Jin, Science 285, 1703 (1999). 

[2] J. D. Weinstein et al, Nature 395, 148 (1998). 

[3] H. L. Bethlem et al, Nature 406, 491 (2000). 

[4] J. Weiner, V. S. Bagnato, S. Zilio, and P. S. Julienne, 

Rev. Mod. Phys. 71, 1 (1999). 
[5] U. Schloder, C. Silber, T. Deuscble, and C. Zimmermann, 

Phys. Rev. A 66, 061403 (2002). 
[6] R. Wynar et al., Science 287, 1016 (2000). 
[7] J. M. Gerton, B. J. Frew, and R. G. Hulet, Phys. Rev. A 

64, 053410 (2001). 
[8] C. McKenzie et al, Phys. Rev. Lett. 88, 120403 (2002). 
[9] M. Pichler and R. Hulet, (private communication). 



3 



[10] E. Timmermans, P. Tommasini, M. Hussein, and A. Ker- 
man, Phys. Rep. 315, 199 (1999). 

[11] S. L. Cornish et al, Phys. Rev. Lett. 85, 1795 (2000). 

[12] T. Weber et al., Science 299, 292 (2002). 

[13] E. A. Donley et al, Nature (London) 417, 529 (2002). 

[14] A. G. Truscott et al, Science 291, 2570 (2001). 

[15] F. Schreck et al, Phys. Rev. A 64, 011402 (R) (2001). 

[16] G. Modugno et al, Science 287, 2240 (2002). 

[17] For a review of the role of Bose enhancement in photoas- 
sociation, see M. Mackie and J. Javanainen, J. Mod. Opt. 
47, 2645 (2000). 

[18] For a brief review of cooperative fermion behavior in four- 
wave mixing, see J. Javanainen, Nature 412, 689 (2002). 

[19] J. Javanainen and M. Mackie, Phys. Rev. A 59, R3186 
(1999). 

[20] D. J. Heinzen, R. Wynar, P. D. Drummond, and K. V. 
Kheruntsyan, Phys. Rev. Lett. 84, 5029 (2000). 

[21] A. Vardi, V. A. Yurovsky, and J. R. Anglin, Phys. Rev. 
A 64, 063611 (2001). 

[22] F. A. van Abeelen and B. J. Verhaar, Phys. Rev. Lett. 
83, 1550 (1999). 

[23] E. Timmermans, P. Tommasini, R. Cote, M. Hussein, 
and A. Kerman, Phys. Rev. Lett. 83, 2691 (1999). 

[24] V. A. Yurovsky, A. Ben-Reuven, P. S. Julienne, and C. 
J. Williams, Phys. Rev. A 60, R765 (1999). 

[25] L. Allen and J. H. Eberly, Optical Resonance and Two- 
Level Atoms, (Dover, New York, 1987). 

[26] S. J. J. M. F. Kokkelmans and M. J. Holland, Phys. Rev. 
Lett. 89, 180401 (2002). 

[27] M. Mackie, K.-A. Suominen, and J. Javanainen Phys. 
Rev. Lett. 89, 180403 (2002). 

[28] T. Kohler, T. Gasenzer, and K. Burnett, Phys. Rev. A 
67, 013601 (2003). 

[29] V. A. Yurovsky and A. Ben-Reuven, Phys. Rev. A 67, 
043611 (2003); M. Mackie et al, (cond-mat/0210131). 

[30] A two-color photoassociation scheme, whereby electron- 
ically excited molecules are transferred to a stable 
molecule state, could of course avert such losses, and 
the work herein indicates that "Bose-stimulated" Raman 
adiabatic passage [M. Mackie et al, Phys. Rev. Lett. 84, 
3803 (2000)] will in principle create for Fermi molecules 
[M. Mackie et al, (physics/0305057)]. 

[31] A. Simoni et al, Phys. Rev. Lett. 90, 163202 (2003). 

[32] E. Timmermans, Phys. Rev. Lett. 87, 240403 (2001). 

[33] M. G. Moore and P. Meystre, Phys. Rev. Lett. 86, 4199 
(2001). 

[34] W. Ketterle and S. Inouye, Phys. Rev. Lett. 86, 4203 
(2001). 



c 

.2 0.8 



£0.6 

(5 

3 0.4 

O 

0) 

O 0.2 



/ \ 






/ \ 


■ (a) / \ 


/ 




N=1 

N=5 




i 


l \\ / 




N=10 




1 

i 


\ V X / 




\ 


1 

- I 


\ X \ / 

\ V 


/A J 
' '" \ / 










V/ 



10 

fit 




FIG. 1. Rabi-like oscillations in the fraction of Fermi-Bosc 
atoms converted to Fermi molecules. The initial number of 
bosonic and fermionic atoms are equal, i.e., Nb = Nf = N. 
(a) The oscillations are complete for an initial number of 
atoms N — 1, while for N = 5 many-body effects lead to fre- 
quency-shifted oscillations that are incomplete and collapse. 
Better short-time agreement with the N = 1 result is obtained 
for N = 10. (b) The oscillations eventually revive. 
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FIG. 2. Rapid adiabatic passage from a Fermi-Bose gas 
of atoms to Fermi molecules. The detuning is swept as 
S(t) = -£f2 2 t, and Nb = Nf = N. (a) For borderline adi- 
abaticity, £ = 1, increasing the number of initial atoms from 
N = 1 to N = 5 indicates that many-body effects reduce the 
efficiency, (b) For N — 5 and £ = 0.1, near-unit conversion is 
still possible, despite many-body effects. 
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